Higher Landau levels contribution 
to the energy of interacting electrons in a quantum dot 
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Properly regularized second-order degenerate perturbation theory is applied to compute the con- 
tribution of higher Landau levels to the low-energy spectrum of interacting electrons in a disk-shaped 
quantum dot. At "filling factor" near 1/2, this contribution proves to be larger than energy differ- 
ences between states with different spin polarizations. After checking convergence of the method in 
small systems, we show results for a 12-electron quantum dot, a system which is hardly tractable 
by means of exact diagonalization techniques. 
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I. INTRODUCTION 

The relevant energy scales entering the Hamiltonian 
of an TV-electron system in a quantum dot (qdot) and 
a magnetic field are the following: the cyclotronic en- 
ergy hw c ~ B, the dot confinement energy huj , and 
the Coulomb characteristic energy e 2 /(/cZ.e) ~ \f~B. In 
strong enough fields, the spacing between Landau lev- 
els (LLs), given by fkv c , is much greater than any other 
scale, and one can restrict the Hilbert space to functions 
built on one-particle states from the first LL. This is the 
ILL approximation 1 , which has been widely used to ob- 
tain exact solutions 2 , to construct the famous v = 1/3 
FQHE functions^, later extended to other filling factors 
by means of the Composite Fermion recipe^ and, in gen- 
eral, has been used to numerically diagonalize the inter- 
acting Hamiltoniani. 

The inclusion of higher LLs in numerical calculations 
turns out to be prohibitive, even for relatively "small" 
systems. Consider, for example, N=12 electrons in a 
qdot at "filling factor" near 1/2, i.e. when the angular 
momentum of the electron droplet is L = —132. Out of 
only 78 one-particle states (orbitals) in the ILL, one can 
construct 674585 Slater determinants, which conform the 
truncated basis for the 12-particle system in the ILL ap- 
proximation. Taking 78 additional orbitals from each of 
the next two LLs causes the basis dimension to be raised 
to more than 172 millions, and the diagonalization of the 
Hamiltonian matrix becomes a very hard computational 
task. 

In the present paper, we show that a way to circum- 
vent the diagonalization of these large matrices is the use 
of properly renormalized degenerate perturbation the- 
ory (PT). We stress that, unlike Monte Carlo and other 
methods focusing on the properties of a particular state, 
by means of PT we obtain, in a single calculation, an ap- 
proximation to the energy spectrum and the correspond- 
ing wave functions of the system. 



The interest in computing the higher LL contribution 
to the energies relies in the fact that, for intermedi- 
ate filling factors, this contribution may be even larger 
than energy differences between states with different spin 
polarizations^. Thus, a correct description of spin exci- 
tations in a system of interacting electrons should take 
account of higher LL effects. Recent work on the issue 
of spin excitations in qdots& has stressed the importance 
of the second LL at v w 2, but at lower v the higher LL 
effects are commonly ignored. 

The plan of the paper is as follows. In the next two 
sections a brief summary of PT and its regularization 
by means of Shank extrapolantsi and the Principle of 
Minimal Sensitivity^ is included for completeness. For 
simplicity, only spin-polarized systems will be studied, 
but any other spin-polarization sector may be treated as 
well. Section llVl is devoted to the results. The 2- and 6- 
electron dots are used as benchmarks where regularized 
second-order PT (PT2) is compared with exact or vari- 
ational results. After validation, the method is applied 
to the 12-electron system mentioned above. Concluding 
remarks are given at the end of the paper. 



II. DEGENERATE PERTURBATION THEORY 

The ILL approximation can be seen, from another 
point of view, as first-order degenerate perturbation the- 
ory. In fact, writing the Hamiltonian in the form: H = 
Ho + V, where H describes free (spin-polarized) elec- 
trons in a magnetic field, and V = V con f + V cou i accounts 
for the external confinement and Coulomb interactions, 
the Hamiltonian matrix in the ILL approximation, 



H 



(i) _ 



(Si\H\Si) = Eobj + (Si\V\S. 



31 ' 



(1) 



where Eq = Nhw c /2 and Si, Sj are Slater determinants 
made up from ILL orbitals, may be seen as the secu- 
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lar matrix of first-order degenerate perturbation theory 9 . 
The degeneracy subspace is spanned by the Si. 

Corrections to are computed in the standard form 9 -. 
The second-order matrix is given by: 

where the sum runs over eigenfunctions of Hq in the or- 
thogonal subspace, {Z\S % ) = 0, and E (Z) = (Z\H \Z). 

We will use Eq. (J2J) to compute higher LLs contri- 
butions to the energy spectrum of an N-electron qdot. 
Notice that the dimension of the secular equation is 
not increased by the inclusion of the second-order cor- 
rections. For the largest systems, an energy cutoff, 
Eq(Z) — Eq < K cut ftiu c , will be imposed to limit the 
number of states entering the sum. We will show results 
with K cut = 2, i. e. three LLs will be included. 

A. The orthogonal subspace 

One can explicitly use the fact that V con f and V cou i 
are, respectively, one- and two-body operators, and ex- 
ploit their symmetries (conservation of total angular mo- 
mentum) in order to carry out the sum only over inter- 
mediate states, Z, having nonvanishing matrix elements 
with one of the external Slater functions, for example 
(Z\V\Sj) ? 0. 
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to higher LLs. And, finally, functions Z% in which one 
occupied level of Sj is moved to an empty level in the ILL 
and a second occupied orbital is moved to a higher LL 
shall also be included. In the first case, both the matrix 
elements of V con f and V cou i could be nonzero, whereas 
in the later two cases only V cou i could have nonvanishing 
matrix elements. 



III. REGULARIZATION OF THE 
PERTURB ATIVE SERIES 

To renormalize the perturbative series (usually an 
asymptotic series) many recipes have been invented. In 
the present paper, we will try Shank extrapolants^ and 
the principle of minimal sensitivity (PMS)£. A variant of 
the later procedure has been recently applied to compute 
the correlation energy of the Coulomb gasiS. 

A. Shank extrapolants 

Shank extrapolants^ are designed to accelerate the con- 
vergence of numerical series. For any three contiguous 
values, Ei, E i+1 and E i+2 , we define the extrapolant: 

P i _ jjffi+g ~ E f+1 ,gN 

Ei + E i+2 — 2E i+ i 

From the series of extrapolants (which will be called 
first order) one can construct the second-order extrap- 
olants, etc. In our case, we have only three values of 
energy, Eq, E\ and E 2 , obtained from PTO, PT1 and 
PT2, and thus there will be only one extrapolant, Fq. 

B. The principle of minimal sensitivity 




FIG. 1: Slater functions in the degeneracy subspace, Sj, and 
in the orthogonal subspace, Z, entering the sum in Eq. 
The example is for a 4-electron system. 

In Fig. ^ we have illustrated this statement for the 
simple 4-electron system. The top of the figure shows 
the occupation corresponding to a given Sj. Then, the 
sum will contain functions Z\ , where one occupied orbital 
of Sj is raised to an orbital in a higher LL (with the 
same angular momentum, I). The sum will also contain 
functions Z 2 , where two occupied levels of Sj are moved 



The PMS starts from the obvious fact that if the 
Hamiltonian is made to depend on an artificial trivial 
parameter a: H — > H(a), then the exact eigenvalues will 
satisfy the equations: 



The perturbative expansion may not, however, respect 
these constraints. The PMS states that an optimal choice 
for a at a given perturbative order is the value at which 
Eq. Q is satisfied. 

In our iV-electron system, described by the Hamilto- 
nian: 
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the introduction of trivial terms and a scaling of coordi- 
nates lead to: 



H = 
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~~2 2/ + ~V _ 2fi~^^~ 



(6) 



where Sw^, is an artificial parameter (magnetic field) and 
L = J2i h- The first sum in Eq. © will be taken as the 
unperturbed Hamiltonian, Hq, and the subsequent terms 
as the perturbation, V. Notice that the term (hw c /2 — 
hui' c /2)L in V will give no contribution to the second- 
order correction because [Ho, L] = 0. 

The parameter Hlu' c will be fixed from the condition 
(0J. Working whithin PT1, one can interpret Hui' c as a 
variational parameter. In PT2, however, the point in 
which Q is satisfied in our examples is a local maximum, 
as will be seen below. 



IV. RESULTS 

The qdot parameters to be used are borrowed from 
Ref. 0. GaAs effective mass, m = 0.067mo, and relative 
dielectric constant k = 12.5 are employed. The confine- 
ment potential is parabolic with Hluq — 3 meV. The bare 
Coulomb interaction is weakened by a factor 0.8 to ap- 
proximately account for quasi-bidimensionality (instead 
of exact bidimcnsionality). We will only deal with spin- 
polarized systems, thus the Zeeman energy is a trivial 
constant and will not be included. 



A. Two electrons 

In the two-electron dot, calculations may be carried 
out semi-analytically. The results for the ground-state 
energy of the triplet state with L = — 1 at £? = 8 T are 
the following: E = 13.82, E 1 = 20.28, and E 2 = 20.16 
meV. The exact energy is E — 20.17 meV. Notice that 
the higher LL contribution to E is -0.12 meV, and that 
E2 is very close to the exact value. 

Regularization of the perturbative series by means of 
the Shank extrapolant yields a value Fq which practi- 
cally coincides with E2- The difference lies in the fifth 
significant figure, not written above. 

On the other hand, as a function of the PMS param- 
eter, B' , we obtain the curves shown in Fig. [21 for the 
magnitudes E X {B') and E 2 (B'). 

Note that the position of the minimum in E\ (B r ) prac- 
tically coincides with the position of the maximum in 
E2{B'). These are the physically acceptable values of 
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FIG. 2: Ei and E2 as a function of B' for the two-electron 
system at B = 8 F. 



B' according to Q). Note also that the PMS-regularized 
ground state energy has the same significant figures as 
the unregularized second-order result, E2(B). 



B. Six electrons 

The exact diagonalization (variational) results of Ref. 
H show that at B = 8 T (fill ing factor near 1/2) the 
contribution of the second and third LLs to the energy 
eigenvalues is around -0.4 meV, i. e. around -0.06 meV 
per electron, as in the two electron system. This magni- 
tude is greater than the energy differences between states 
with different spin polarizations. 

The lowest spin-polarized states in each angular mo- 
mentum tower (the yrast spectrum) are shown in Fig. 
[3] Energy jumps between adjacent angular momentum 
states are about 0.6 meV. The PT2 results are shown 
to lay 0.04 meV below the variational ones, that is, al- 
most inside the estimated error bars which are about 0.03 
meV. These PT2 results are obtained with K cut — 2, as 
mentioned above. 

Regularization of the perturbative series leads to re- 
sults very similar to those in the two-electron dot: the 
significant figures (in this case five) are not changed. This 
means that one can in practice estimate the ground state 
energy from the PT2 result. 



C. Twelve electrons 

We show in Fig. 01thc perturbative low-lying spectrum 
for the 12-electron dot. In this situation, the pertur- 
bative results are the only ones at our disposal because 
exact diagonalization is practically impossible. The mag- 
netic field is 10 T in this figure. Notice the energy jumps 
between adjacent angular momentum states, as in the 6- 
electron dot. The contribution of the higher LLs to the 
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FIG. 3: Variational and PT2 results for the low-lying states 
in the 6-electron dot. B — 8 T. 
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FIG. 4: PT1 and PT2 estimates of the low-lying spin- 
polarized energy levels of twelve electrons. B — 10 T. 



are reliable as well. The generation of the second-order 
matrix, Eq. for a single L value takes around one 
week CPU time in a small computer cluster with ten 
processors at 2.4 GHz. The PT2 matrix is a factor of 
ten less sparse than the PT1 matrix, occupying around 
15 GB of hard disk. Consequently, diagonalization by 
means of a Lanczos algorithm takes a factor of ten more 
time. 



V. CONCLUDING REMARKS 

We have shown the feasibility of computing the energy 
eigenvalues and eigenfunctions of relatively large elec- 
tronic quantum dots in a magnetic field by using second- 
order degenerate perturbation theory to take account of 
the contribution of the higher LLs. This contribution 
proves to be of the same order of energy differences be- 
tween states with different spin polarizations. The con- 
structed effective Hamiltonian in the ILL can be diag- 
onalized by means of standard Lanczos techniques. We 
presented calculations for the spin-polarized states of a 
12-electron quantum dot at filling factor near 1/2. When 
three LLs are included, the dimension of the truncated 
Hilbert space for this problem is larger than 170 millions. 

From a more general point of view, a similar scheme 
could be applied to any problem with well differentiated 
energy scales, in which one could identify, at least ap- 
proximately, degeneracy subspaces. Examples of such 
problems are the study of mixing between vibrational 
and electronic degrees of freedom in relatively large 
molecules, or the study of mixing between intra- and 
inter-shell excitations in relatively large nuclei. 
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